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ABSTRACT 

This  document  describes  a  procedure  for  designing  recursive 
digital  filters  from  continuous-time  filters  when  the  ratio  of  the 
sampling  frequency  to  the  pole  frequency  is  small.  The  discrete-time 
transfer  function  of  digital  filters  obtained  by  this  method,  called 
the  acceleration-invariant  transformation,  has  been  derived  from  the 
partial  fraction  expansion  of  the  continuous-time  transfer  function. 

An  error  function  between  the  digital  and  analog  transfer  functions  has 
also  been  deduced.  Finally,  the  performance  of  the  proposed  method  is 
demonstrated  by  plotting  the  frequency  response  of  high-order 
Butterworth  and  elliptic  filters./ _ 

RESUME 


Ce  document  dficrit  une  m§thode  pour  concevoir  des  filtres  num6- 
riques  rScursifs  A  partir  des  filtres  analogiques  lorsque  le  rapport 
entre  le  taux  d'dchantillonnage  et  la  frequence  du  pdle  est  faible.  La 
fonction  de  transfert  en  z  des  filtres  num£rlques  obtenus  par  cette 
mSthode,  appelAe  1* invariance  A  1' acceleration,  a  6te  dedulte  de  la 
decomposition  en  fraction  partielle  de  la  fonction  de  transfert  des 
filtres  analogiques.  Une  fonction  qui  permet  de  calculer  l'erreur 
entre  la  fonction  de  transfert  numerique  et  la  fonction  de  transfert 
analoglque  a  aussi  6te  dedulte.  Finalement,  la  performance  de  cette 
methode  est  demontree  en  traqant  la  reponse  en  frequence  des  filtres 
Butterworth  et  elliptiques  A  ordres  eiev6s. 
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1.0  INTRODUCTION 


This  document  describes  an  alternative  approximation  method  for 
recursive  digital  filters  when  the  ratio  of  the  sampling  frequency  to 
the  cutoff  or  center  frequency  (fg/fc)  is  small.  Recursive  digital 
filters  are  usually  designed  by  using  the  impulse-invariant  method  and 
the  bilinear  transformation.  These  methods,  which  constitute  textbook 
material  (Refs.  1  to  5),  permit  the  derivation  of  a  digital  filter  from 
a  suitable  analog  filter  satisfying  a  set  of  prescribed  specifications. 
They  also  provide  valid  matches  with  their  corresponding  analog  filters 
when  the  ratio  1*  sufficiently  high.  At  low  f^/f^,  these  methods 

have  serious  drawbacks  that  lead  to  inappropriate  approximations. 

The  standard-z  transformation  (also  called  the  impulse-invariant 

transformation)  is  a  method  whereby  the  response  of  the  derived  digital 

filter  to  an  impulse  is  identical  to  the  sampled  Impulse  response  of 

the  continuous-time  filter.  This  method  yields  valid  matches  with  the 

corresponding  analog  filter  only  at  high  sampling  rates  and  it  is 

satisfactory  when  the  analog  system  is  sufficiently  band-limited  prior 

to  transformation.  At  low  ratios  of  f  /f  ,  this  method  can  be  applied 

sc 

only  to  transfer  functions  in  which  the  denominator  degree  exceeds  that 
of  the  numerator  by  at  least  two.  The  aliasing  effects  obtained  by 
this  transformation  render  it  entirely  useless  for  digital  high-pass 
and  band-reject  filters.  In  addition,  for  elliptic  low-pass  and 
band-pass  filters,  the  equiripple  character  of  the  stop  band  is 
destroyed  by  the  aliasing  effect.  The  bilinearly  transformed  filters 
are  essentially  identical  to  the  original  analog  filter  when  the  ratio 
f8/fc  is  greater  than  10.  However,  this  transf ormation  results  in  a 
nonlinear  relationship  between  the  analog  and  the  digital  frequencies. 
It  implies  that  the  frequency  response  of  the  continuous  system  being 
transformed  will  be  a  warped  version  of  the  analog  frequency  response. 
This  effect,  which  is  observed  particularly  when  the  ratio  fg/fc  Is 
smaller  than  10,  renders  the  method  quite  complicated  when  designing 
adaptive  filters.  In  addition,  neither  the  impulse  response  nor  the 
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phase  response  of  Che  analog  filter  Is  preserved  In  a  digital  filter 
obtained  by  the  bilinear  transformation. 

Other  methods  of  deriving  recursive  forms  of  continuous-time 

transfer  functions  of  analog  filters  are  an  extension  of  the  concept  of 

the  impulse  invariance.  In  the  sinusoid  invariance  (Ref.  6),  the  step 

invariance  (Refs.  7  and  8)  and  the  ramp  invariance  (Ref.  8),  the 

digital  filter  design  is  optimized  for  a  sequence  of  sinusoids,  steps 

or  ramps  joining  the  sampled  values  as  an  input  signal.  These  methods 

attenuate  the  aliasing  effect  encountered  in  the  impulse  invariance. 

The  sinusoid-invariant  method  preserves  the  phase  characteristics  of 

analog  filters  and  it  is  suitable  to  the  design  of  filters  with  a 

constant  group  delay.  The  step-invariant  method  can  be  applied  to 

transfer  functions  in  which  the  degree  of  the  denominator  can  exceed 

that  of  the  numerator  by  at  least  one.  In  the  ramp-invariant  method, 

the  numerator  degree  can  be  as  high  as  the  denominator  degree.  In 

Ref.  8,  it  is  shown  that  the  step  and  ramp  invariances  of  low-pass  and 

band-pass  filters,  if  they  are  realized  in  a  parallel  form,  produce 

digital  filters  whose  frequency  response  does  not  depend  on  the 

sampling  frequency  in  the  frequency  band  of  0  to  f g/2.  The  frequency 

response  of  the  step-invariant  high-pass  and  band-reject  filters 

realized  in  a  parallel  or  a  cascade  form  is  unacceptable  for  small 

ratios  because  of  the  aliasing  effects.  The  ramp  invariance 

realized  in  a  cascade  form  can  be  used  to  digitize  high-pass  and 

band-reject  filters.  The  ratio  f  If  can  be  decreased  below  10  for 

sc 

Butterworth  filters  but  the  equiripple  character  in  the  stop-band 
response  of  elliptic  filters  is  preserved  by  a  ratio  of  at  least  10. 

The  approximation  method  described  in  this  report  is  derived 
from  analog  filters  such  that  an  acceleration  invariant  response  is 
maintained.  The  discrete-time  transfer  function  of  the  digital  filters 
is  obtained  from  the  partial  fraction  expansion  of  the  continuous-time 
transfer  function  to  lead  directly  to  a  realization  in  a  parallel  form. 
An  error  function  between  the  transfer  functions  of  the  digital  and 
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analog  filters  has  also  been  deduced.  This  method  is  shown  to 

approximate  correctly,  on  a  wide  bandwidth,  the  frequency  responses  of 

low-pass,  band-pass,  high-pass  and  band-reject  analog  filters  for  small 

ratios  of  f  If  . 

s  c 

This  work  was  performed  at  DREV  between  November  1982  and  April 
1983  under  PCN  21J05,  Guidance  and  Control  Concepts. 

2.0  DERIVATION  OF  THE  ACCELERATION- INVARIANT  METHOD 


Assume  that  the  coefficients  of  an  analog  filter  with  Laplace 
transformation  (i.e.  transfer  function)  given  by 


V> 


I 


«•>  f ./ 

i«o 


[1] 


are  known.  Y(s)  and  X(s)  are  respectively  the  Laplace  transformation 
of  the  output  y(t)  and  the  input  x(t).  M'  must  be  smaller  or  equal  to 
N'.  If  we  assume  that  all  real  and  complex  poles  are  simple,  eq.  1  can 
be  rewritten  in  its  partial  fraction  expansion  as 


H(s) 


A  + 
o 


N 

l 

i-1 


[2] 


Equation  2  covers  low-pass,  high-pass,  band-pass  and  band-reject 
filters  among  the  well-known  classes  such  as  the  Butterworth,  Bessel, 
Chebyshev  and  elliptic  filters. 


The  s  domain  filter  response  of  eq.  2  to  an  acceleration  input 
is  written  as 


ya(b)  "  V8)x(s)  “  V8)/®3 


[3] 
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Substituting  eq.  2  into  eq.  3  and  finding  the  inverse  Laplace 
transformation  of  this  relation  yield  the  following  time  function: 


ya  (t) 


{V!  -  v  +  M1-  V))  [4] 


2  i-1  2Pl  P2 


Hence,  the  sampled  sequence  of  YA(t)  can  be  obtained  by  applying 
the  z  transformation  to  eq.  4.  This  gives: 

r2„-l,i  j.  N  A ,T2z-1(1  +  z“l) 


A0T^z-‘(l  +  z“i) 

V->-  V(V  - 


A1Tz-1 


I  ^ 


i-1  2Pi(l  -  z~l)3 

-P.T 

A1z‘1(l  ‘  e  ) 


P2(l  -  z“l)3  p3(l  -  z_1)(l  -  z-1e  PiT) 


[5] 


Ts 

where  z=  e  and  is  the  sampling  period.  Yp(z)  can  also  be  expressed 
as 


Vz)  "  hd(z)-x(z) 


[6] 


where  HD(z)  is  the  discrete-time  transfer  function  of  the  digital 
filter.  X(z),  the  sampled  sequence  of  the  acceleration  input,  is 


X(z) 


T^-^l  +  z~l) 
2(1  -  z-1)3 


[7] 


Hence,  H_(z)  is  derived  by  combining  eqs.  5,  6  and  7.  This 
operation  yields 


»D<Z> 


lisl 

X(z) 


.  N  C. z“2  +  D,z_1  +  E, 

A  ,  r  1  i  i 

0  +  i  - iTt - 

i-1  (1  +  z_1)(l  -  e  1  z_I) 


[8] 
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A  2A  -P  I  2A  -P  T 

where  C.  -  -  +  -±)  e  +  (1  +  e  ) 

1  P  P^T  P3T^ 

i  i  1 


-PI 

A.  4A  -P  T  2A .  (1  +  e  ) 

Di  -  <r  -  (1 " e  > +  — 

i  i  i 


_  _  i  i  .  1  ,,  i  . 

Ei  P  p2t  +  p3T2  ~  6  5 

A  digital  filter  whose  response  to  an  acceleration  excitation 
corresponds  to  that  of  the  analog  filter  has  then  been  obtained. 

3.0  ERROR  FUNCTION  BETWEEN  THE  ANALOG  AND  THE  ACCELERATION- 
INVARIANT  DIGITAL  FILTER 


The  purpose  of  this  chapter  is  to  establish  an  error  function 
between  the  continuous -time  and  the  discrete-time  transfer  functions. 

A  relation  can  be  found  by  regarding  the  response  of  the  digital  filter 
as  the  sampled  response  of  the  corresponding  analog  filter.  The 
Fourier  transform  of  y(t)  is  obtained  from  Ref.  1  as 


V 


>»T) 


where  u>g  »  2x/T 


The  substitution  of  eq.  9  into  eq.  3  gives 


v«>T> 


1  7  ha  (Jm  *  J-V 

*  m—  +  “s> 


[9] 


[10] 
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The  response  of  the  digital  filter  to  an  acceleration  input 
excitation,  which  can  be  derived  from  the  relation 
(Vz>  “  »D(z).X(z)),  gives 


>T  HD(e>T)T2e-^T(l+e-^T) 


aci-e"^1)3 


[11] 


A  relation  between  the  analog  and  digital  transfer  functions  is 
obtained  by  substituting  eq.  11  into  eq.  10.  This  yields 


HD(ej“T)_  _  2(1  -  e"*1*1) 

jT3e-J“T(l  +  e'J011)  (u  +  wflg)3 


-  j-y  2] 


For  sufficiently  large  «*»  _,  the  effect  of  frequency  aliasing  can 
be  neglected.  Under  this  condition,  eq.  12  can  be  written  as 


-  E(eJ“T>HA(>) 


[13] 


IwT 

where  E(eJ  )  is  an  error  function  given  by 


— -fcii  T  ^ 

E^1) _ -■£■■■> 

v  '  ,  ,  -*uT 


jw3T3e'JwT(i  +  e  JuT) 


[14] 


The  magnitude  and  phase  responses  of  the  error  function  (eq.  14) 

have  been  plotted  in  Fig.  1  from  fg/101*  to  fg/2.5.  It  is  shown  that 

the  magnitude  response  of  the  digital  filter  is  Identical  to  that  of 

its  analog  counterpart  over  the  frequency  band  of  fg/104  to  fg/10.  The 

error  in  magnitude  Increases  gradually  to  attain  2.5  dB  at  fg/2.5.  No 

error  in  phase  is  observed  from  f  /104  to  f  /2.5.  From  f  /2.5,  the 

8  8  s 

error  increases  to  attain  *  at  f  /2.  The  acceleration-invariant 

8 

filters  are  thus  limited  to  the  frequency  range  of  0  to  fg/2.5. 
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Magnitude (dB)  Phase(deg) 


FIGURE  1  -  Frequency  response  of  the  error  function  between  the  analog 
and  digital  transfer  functions 


4.0  PERFORMANCE  OF  THE  ACCELERATION-INVARIANT  METHOD 
FOR  HIGH-ORDER  BUTTERWORTH  AND  ELLIPTIC  FILTERS 


Figures  2  to  5  show  the  magnitude  and  phase-frequency  responses 
of  high-order  Butterworth  and  elliptic  digital  filters  obtained  by  the 
acceleration-invariant  method.  In  all  examples,  the  response  of  digit¬ 
al  filters  is  presented  for  various  sampling  frequencies  (100,  10  and 
4  kHz)  and  it  is  accompanied  by  the  response  of  the  corresponding  ana¬ 
log  filters.  Figures  2  and  3  illustrate  fifth-order  low-pass  and  high- 
pass  Butterworth  and  elliptic  filters.  The  cutoff  frequency  (fc)  and 
the  gain  of  these  filters  are  maintained  at  1  kHz  and  1  respectively. 

In  addition,  for  elliptic  filters,  the  stop-band  attenuation  is  at 
least  40  dB  and  the  pass-band  ripple  is  limited  to  3  dB.  The  band-pass 
and  the  band-reject  Butterworth  and  elliptic  filters  of  Figs.  4  and  5 


I 
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were  derived  from  a  transformation  performed  on  the  corresponding 
low-pass  filters  (Ref.  1).  The  lower  and  upper  cutoff  frequencies  of 
these  filters  are  located  at  700  Hz  and  1.0  kHz.  Their  center  frequen¬ 
cy  is  noted  by  fc-  The  poles  and  zeros  of  analog  filters  are  given  in 
Appendix  A. 

The  acceleration  invariance  of  the  low-pass  and  band-pass 

filters  produces  digital  filters  whose  frequency  response  matches  very 

closely  that  of  their  analog  counterpart  in  the  frequency  band  of  0  to 

fs/2.5.  In  this  region,  the  frequency  response  does  not  depend  on  the 

sampling  frequency.  For  elliptic  filters,  the  equiripple  character  of 

the  stop-band  response  is  correctly  preserved.  The  high-pass  digital 

filters  are  essentially  identical  to  the  original  analog  filter  when 

the  ratio  f^/^  is  greater  than  or  equal  to  10.  When  it  decreases 

below  10,  the  aliasing  effect  can  be  observed  in  the  magnitude  response 

of  the  digital  high-pass  filters.  The  Butterworth  filter  exhibits  an 

error  in  magnitude  and  in  phase  in  the  low-frequency  band.  However, 

when  f  is  set  to  4f  ,  this  error  does  not  exceed  10  dB  in  magnitude 

and  20°  in  phase  at  f^/10.  The  high-pass  elliptic  filter  does  not 

perfectly  preserve  the  equiripple  character  of  the  stop-band  response 

when  i^/f^  is  smaller  than  10.  The  phase  response  is  also  slightly 

disturbed  around  the  resonant  frequency.  The  band-reject  digital 

filters  react  in  the  same  manner  as  the  high-pass  filters.  Their 

frequency  responses  coincide  with  those  of  the  corresponding  analog 

filters  when  the  ratio  f  /f  is  greater  than  or  equal  to  10.  When  it 

sc 

becomes  smaller  than  10,  the  attenuation  of  the  Butterworth  filter, 

the  phase  response  and  the  equiripple  character  of  the  elliptic  filter 

in  the  stop-band  are  not  correctly  reproduced.  Although  all  the 

characteristics  of  the  digital  high-pass  and  band-reject  filters  are 

not  preserved  when  the  ratio  f  /f  is  set  to  a  value  smaller  than  10, 

sc 

no  gain  is  Introduced  by  the  acceleration-invariant  transformation. 
Furthermore,  this  method  does  not  move  the  poles  of  the  digital  filters 
as  a  function  of  the  sampling  frequency  and  it  can  be  applied  to  a  wide 
frequency  band. 


I 
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FIGURE  3  -  Frequency  responses  of  Butterworth  and  elliptic  band-pass 
digital  filters  obtained  by  the  acceleration-invariant 
method  for  various  ratios  f8/fc  (4,  10,  100) 
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FIGURE  A  -  Frequency  responses  of  Buttervorth  and  elliptic  high-pass 
digital  filters  obtained  by  the  acceleration-invariant 
method  for  various  ratios  fg/fc  (A,  10,  100) 
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Frequency  responses  of  Butterworth  end  elliptic  bend-reject 
digital  filters  obtained  by  the  acceleration-invariant 
aethod  for  various  ratios  f8/fc  (4,  10,  100) 
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5.0  CONCLUSION 


An  approximation  method  has  been  described  whereby  a  digital 
filter  can  be  derived  from  an  analog  filter  while  maintaining  an 
acceleration-invariant  response.  The  discrete-time  transfer  function 
of  digital  filters  obtained  by  this  method  has  been  deduced  from  the 
partial  fraction  expansion  of  the  continuous-time  filter  yielding  a 
parallel  realization. 

The  performance  of  the  proposed  method  has  been  determined  by 
giving  the  frequency  response  of  high-order  Butterworth  and  elliptic 
filters.  The  acceleration  Invariance  of  low-pass  and  band-pass  filters 
produces  digital  filters  whose  frequency  responses  match  perfectly 
those  of  their  analog  counterparts.  Also,  these  responses  do  not 
depend  on  the  ratio  f8/fc  *n  the  frequency  band  of  0  to  fs/2.5.  An 
accurate  frequency  response  of  high-pass  and  band-reject  filters  can  be 
obtained  when  the  ratio  f8/fc  is  maintained  equal  to  or  greater  than 
10.  When  it  decreases  below  10,  the  frequency  response  of  the  digital 
filters  remains  acceptable  but  the  equiripple  character  of  the  elliptic 
filter  in  the  stop-band,  the  magnitude  response  of  the  Butterworth 
high-pass  filter  in  the  low-frequency  band  and  the  peak  stop-band 
attenuation  of  the  Butterworth  band-reject  filter  are  not  perfectly 
reproduced  by  the  transformation. 

The  acceleration-invariant  method  is  then  less  sensitive  to  the 
frequency  aliasing  in  comparison  with  the  impulse-invariant,  the  step- 
invariant  and  the  ramp-invariant  methods.  In  addition,  this  method 
exhibits  no  frequency  warping,  a  drawback  met  in  the  bilinear  transfor¬ 
mation.  However,  this  improvement  produces  an  increase  in  the  order  of 
the  acceleration-invariant  filter.  The  numerator  and  denominator 
degrees  of  the  acceleration-invariant  filter  exceed  by  one  those  of  the 
ramp-invariant  or  the  bllinearly  transformed  filter. 
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APPENDIX  A 


Continuous-Time  Transfer  Functions  of 
the  Butterworth  and  Elliptic  Filters 

This  appendix  gives  the  continuous-tiae  transfer  functions  of 
the  Butterworth  and  elliptic  filters  used  in  Chapter  4. 

The  transfer  functions  of  these  filters  are  written  as 

[A.l 

1)  Fifth-order  Butterworth  low-pass  filter: 


H(s)  -  K 


i5l  <8  +  V 


i5i  <*  +  v 


-  Cutoff  frequency  -  1000  Hz 

-  M*  *  number  of  zeros  “  0 

-  N'  -  number  of  poles  -  5 

-  K  -  9.79  x  1018 

s 

-  p1>2  -  1941  ±  j  5975 

-  P3.4  -  5083  ±  j  3693 

-  ps  -  6283 

2)  Butterworth  high-pass  filter  derived  from  the  low-pass 
filter: 

-  Cutoff  frequency  ■  1000  Hz 

-  M'  -  5,  N*  -  5,  K  •  1 

-  poles  are  identical  to  those  of  the  low-pass  filter 

-  zeros  are  all  located  at  0 


[ 


•y 


1 


» 


. 
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3)  Butterworth  band-pass  filter  derived  f roe  the  low-pass 
filter: 

-  lower  cutoff  frequency  “  700  Hz 

-  upper  cutoff  frequency  “  1000  Hz 

-  M*  -  5,  N'  -  10,  Kb  -  2.37  x  1016 

-  zeros  are  all  located  at  0 

-  Pi. 2  *  340  ±  6221 
p3,„  -  242  ±  4428 
p5,6  -  843  ±  5785 
p7,8  -  682  ±  4677 
Pq » 10  -  942  t  5172 

4)  Butterworth  band-reject  filter  derived  f ron  the  low-pass 
filter: 

-  lower  cutoff  frequency  •  700  Hz 

-  upper  cutoff  frequency  ■  1000  He 

-  M’  -  10,  N*  •  10,  K  -  1 

-  poles  are  identical  to  those  of  the  band-pass  filter 

-  zeros  are  located  at  0  t  j  5257 

5)  Fifth-order  elliptic  low-pass  filter: 

-  cutoff  frequency  ■  1000  Hz 

-  M'  -  4,  N'  -  5,  k#  -  240 

-  Bpj  *  0  t  9309,  «3,4  ■  0  t  j  6865 

-  p1>2  -  775.2  t  j  4367.4,  pj,,,  -  155.5  t  j  5802 
PS  -  1479.1 
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6)  Elliptic  high-pass  filter  derived  from  the  low-pass  filter 

-  cutoff  frequency  ■  1  kHz 

-  M*  -  65,  N'  -  5,  K  -  1 

8 

-  mlt2  -  0  t  j  4240,  a3,H  -  0  ±  j  5750 
m5  “  0, 

-  Pl,2  -  182.3  ±  j  6799,  p3M>  -  1555  ±  j  8763 
p5  -  26689 


7)  Elliptic  band-pass  filter  derived  from  the  low-pass  filter 

-  lower  cutoff  frequency  ”  700  Hz 

-  upper  cutoff  frequency  -  1000  Hz 

-  M*  ■  9,  N'  -  10,  Kb  -  1195.2 

-  m1>2  -  0  t  i  6836,  »JPH  »  0  t  j  4043 
■5.6  -  0  1  j  6386,  B7,8  -  0  t  j  4327 

«9  *  0, 

Pi , 2  -  27,2  t  j  6199,  p3,„  -  222  ±  J  5252 
p5,6  -  130  t  j  5951,  p7,8  -  102  ±  j  4641 
P9 . 10  “  10-2  ±  j  4458 

8)  Elliptic  band-reject  filter  derived  from  the  loir-pass 
filter: 


-  lower  cutoff  frequency  -  700  Hz 

-  upper  cutoff  frequency  _  1000  Hz 

-  M'  -  10,  N'  -  10,  Kb  -  1 

-  mi, 2  -  0  1  j  5257,  m3,q  -  0  t  j  5931 
■5.6  *  0  ±  j  4659,  *7,8  -Of  j  6190 
m9,10  ■  0  t  j  4464 

-  pl>2  -  4003  t  j  3407 

P3,„  -  22  1  j  4335,  p5,6  -  32.6  t  j  6375 
p7,8  -  176.6  i  j  4100,  p,,10  -  290  1  j  6729 
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